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We study the superconducting proximity effect on the helical edge states of time-reversal- 
symmetric fractional topological insulators(FTI). The Cooper pairing of physical electrons results 
in many-particle condensation of the fractionalized excitations on the edge. We find localized zero- 
energy modes emerge at interfaces between superconducting regions and magnetically insulating 
regions, which are responsible for the topological degeneracy of the ground states. By mapping 
the low-energy effective Hamiltonian to the quantum chiral Potts model, we determine the opera- 
tor algebra of the zero modes and show that they exhibit nontrivial braiding properties. We then 
demonstrate that the Josephson current in the junction between superconductors mediated by the 
edge states of the FTl exhibit fractional Josephson effect with period as multiples of An. 
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I. INTRODUCTION 

Topological phases are often characterized by gapless 
boundary excitations which do not arise in the low-energy 
part of local lattice Hamiltonians with the same dimen- 
sion. For example, the integer/fractional quantum Hall 
states support chiral edge excitations, which have only 
"half" the degrees of freedom as particles moving in 
a onc-dimensional lattice. Recently discovered Topo- 
logical Insulators(TI) with Time-Reversal(TR) symme- 
try support gapless helical boundary states in both two 
and three dimension^^^'^. These boundary states can be 
gapped out by adding certain symmetry-breaking pertur- 
bations(e.g. superconducting or magnetic order) which 
in many cases leads to exotic phases. A notable exam- 
ple is the Px + ipy superconducting state created by su- 
perconducting proximity effect on the surface states of 
three-dimensional Tl^' which exhibits non-Abelian Majo- 
rana zero modes in vortices. Similar physics can also be 
realized on the edges of two-dimensional quantum spin 
Hall insulator&iS.^ where Majorana zero modes appear 
at the interfaces between superconducting and magnetic 
gapped regions. The Majorana zero modes exhibit un- 
usual properties such as 47r-period Josephson effeclliSHill 
and non-Abelian statistic^SK^^, which have important 
application in quantum information processing''^^ 

On the other hand, strongly correlated topologi- 
cal phases , such as Fractional Quantum Hall(FQH) 
state^l^'^-^, are usually associated with the fractionaliza- 
tion of quantum numbers. It would be even more in- 
teresting to study the quantum phases originated from 
symmetry breaking in the fractionalized boundary states. 
In this paper we study superconducting proximity effect 
on the edge of two-dimensional fractional topological in- 
sulators (FTI)^-''-^, which can be regarded as the TR 
symmetric generalization of Laughlin states with filling 
fraction v = We focus on the properties of the gap- 
less edge states of FTI brought in contact with an s-wave 
superconductor. 

Our main finding is that the electron fractionalization 
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FIG. 1: Schematic illustration of superconductor-ferromagnet 
junctions on the edge of 2D FTI. Localized zero modes appear 
at the interface between the superconducting regions(SC) and 
the magnetically insulating regions(FM). 



drastically changes the superconducting proximity effect 
on the edge states, as opposed to the non-interacting 
quantum spin Hall insulators. Most remarkably, domain 
walls between the superconducting and magnetic gapped 
regions are found to carry localized zero modes with 
quantum dimension d — yjlra. We also determine the 
operator algebra satisfied by the zero modes which can be 
regarded as a generalization of Majorana fermions. 
We then argue that the zero modes exhibit non-Abelian 
statistics upon adiabatic exchanging and determine the 
braiding matrices. We also discuss possible physical in- 
dications of the unusual zero modes and propose an un- 
conventional Josephson effect with period 47rr7i. 

The paper is organized as follows: in Sec. II we review 
the effective edge theory of S'^-conserved FTI. In Sec. HI 
and IV we study the gapped phase driven by s-wave pair- 
ing. We reveal the topological degeneracy of the ground 
states and find the localized zero modes. In Sec. V we 
discuss the braiding of the zero modes. 
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II. EDGE THEORY OF FTI 

We start by reviewing the effective edge theory of 
FT F3 | 24 | 26 | mainly consider FTI with conserved 
where the two spin species each form Laughhn states with 
fiUing fraction v = — where m is an odd integer, under 
a spin-dependent magnetic field B = Bozaz- It supports 
gapless edge states described by a helical Luttinger liq- 
uid model. Levin and SterrP^ have recently shown that 
such edge states are protected by TR symmetry if and 
only if (tsh/c* is odd. Here ugH is the spin Hall conduc- 
tance measured in units of e/27r and e* is the elementary 
charge. Here the ratio crsH/e* = 1 implying the robust- 
ness of the edge states protected by TR symmetry. The 
effective Lagrangian dens ity go verning the dynamics of 
the edge states is given by^IMl 

= ^ ^ {K„„,dt<l)„d^<f)„. - V„„,d^(l)„d^(l)^,). (1) 

The K matrix is K = mcFz and V is the renormalized 
charge velocity matrix. The chiral bosonic fields 0^ sat- 
isfy Kac-Moody algebra 

[(t)„{x),(j)„'{x')] = (cr^)^^/— sgn(x - a;'). (2) 

m 

To simplify our derivation, we define if = ^i^Rt ~^ 
4>Li),0 = — 4'B.t)- They then satisfy the canon- 

ical commutation relation [ip{x),dx'0{x')] = ii:5{x — x'). 
The Hamiltonian of the edge theory can be expressed as 

H = Jdx^ [mg{d,9)' + img)-\d,ipf] , (3) 

which describes a Luttinger liquid with Luttinger param- 
eter mg. Here g — I ii — 0, and g > 1(< 1) if 
% < 0(> 0). 

In this bosonic theory, the electron density p^r is ex- 
pressed as = i^dx<j>a — :^dx(p- The physical electron 
creation operators are given by ijj^ = y^^^*""^^"^- Here 
a is a regularization factor. Notice that proper Klein fac- 
tors should be included in the expression to ensure the 
correct fermionic commutation relations but it turns out 
that they are not relevant for our discussion below, so we 
omit them to simpify the notations. We also omit the 
spin indices of the bosonic fields since they are locked to 
the chiral indices. We adopt the convention that under 
TR transformation (f>^ ^ (j)^, (j)^ ^ (j)^ — ^ to guarantee 
that V't ~^ — ^ — Consequently, ^ — > — ^ — jm' 
One can also define V'qp ~ e"^^"^ which creates a charge 
e/m quasiparticle(QP) on the edge. 

The proximity effect of an s-wave superconductor can 
be taken into account by adding the following pairing 
term to the edge theory: 

T^sc = Ae^^VLV-L +h.c.. (4) 



Here A is the induced superconducting gap and x is the 
superconducting phase. Using the bosonic representa- 
tion, we can write the pairing term as 

-Hsc = —e'x^im{<l>L-<l>H) _^ h.c. = — cos(2m6l -I- x)- 
27ra 7ra 

(5) 

We first assume A is a weak perturbation and study 
its fate using perturbative Renormalization Group (RG) 
analysis. The flow of the coupling A under RG is given 
by 

Here £ = ln(///o) is the dimensionless flow parameter. 
The pairing term is relevant when g > y- For m > 3, in 
the "non-interacting" case g = 1 the pairing is irrelevant. 
This is a direct consequence of the electron fractionaliza- 
tion. However, we should remark here that although the 
term is perturbatively irrelevant when g < ^ , it can still 
gap out the edge Luttinger liquid if the bare value of A is 
large enough when the perturbative RG analysis breaks 
down. 



III. TOPOLOGICAL DEGENERACY IN THE 
GAPPED PHASE 

We now consider the gapped phase resulting from the 
mass term cos{2ni6 + x)- We choose the gauge x = 0. 
Semi-classically, 6 is pinned to the minima of A cos 2m9 
located at: 

e„ = — + —,71 = 0,1,. ..,2m-l. (7) 
Zm m 

Therefore there are 2m degenerate ground states which 
we denote by |6'„). We notice that when m = 1 it re- 
duces to the well-known to polog ical degeneracy of one- 
dimensional Major ana chairpSHMl Since under TR 9 goes 
to —9, the semiclassical ground states |0„) are not TR in- 
variant when TO > 1. 

However, it is important to realize that the one- 
dimensional system under consideration is already the 
boundary of a 2D system and thus it must be a closed 
manifold without any boundaries. If the whole edge 
is gapped out by the superconducting proximity effect, 
there is a unique ground state fixed by the boundary 
condition. The ground state degeneracy only occurs in 
an open geometrjl^. To effectively create boundaries on 
the edge, we need to introduce a different mass term, e.g. 
TR-breaking perturbations by applying Zeeman field or 
proximity effect to a ferromagnetic insulator to gap out 
some regions of the edge liquid. See Fig. [T] for an illus- 
tration of the setup. To be specific, we consider inducing 
ferromagnetic order in some region of the edge, by adding 
the following mass term to the Hamiltonian: 

"Hfm = AfmV'bV'l + h-c- = — '-^ cos 2ip. (8) 
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Now consider a setup in which the edge is divided into 
2n segments, n of which are gapped out by the super- 
conducting proximity effect and the others by proximity 
to ferromagnets, arranged in an alternating order. From 
our semiclassical analysis, each superconducting segment 
has 2m ground states. Superficially there would be to- 
tally (2TO)"-fold degeneracy. However, due to the global 
conservation law of QP number mod 2m (see below) , the 
degeneracy is actually reduced to (2m)"~^. Since the 
entire edge is gapped away from the interfaces between 
different regions, the ground state degeneracy can only 
arise from zero modes localized at the interfaces. Because 
there are 2n interfaces, each zero mode has a quantum 
dimension d = \j2m. 

To understand the nature of the degenerate ground 
states, we notice that the total charge density on the edge 
is given by /? = ^^'^xt^s. Its commutation relation with 
B is \p(x),Q{^')\ = - x'). Therefore the following 

relation can be derived: 



IV. LOW-ENERGY EFFECTIVE 
HAMILTONIAN AND ZERO-ENERGY 
BOUNDARY MODES 

In the previous section we use semiclassical analysis to 
study the ground state properties of the gapped phase. 
We argue that it is natural to relate the ground state 
degeneracy to localized zero modes on the SC-FM inter- 
faces where the spectra gap has to close. Heuristically, 
the zero modes are transformations between the ground 
states. In this section we put these considerations on 
a firmer ground by deriving an effective Hamiltonian to 
describe the low-energy quantum fluctuations around the 
semiclassical ground states. 

First it will turn out to be convenient to regularize the 
bosonic Hamiltonian on a lattice (after proper rescaling 
of space): 

B = uS^ |"-gcos(6'i-6',+i)-f — cos2m6'i]. (11) 

^—^ '■ 2 TTM 



m 



(9) 




Here Q — J dx p(x) which counts the number of QP's 
residing on the edge. We can then linearly superpose 
\9j) to obtain eigenstates of e^'^''^; 



c^"^|0,), g|n) (10) 



Here uj = . Therefore, the 2m degenerate ground 
states have different mod 2m QP number. This is anal- 
ogous to the two-fold ground state degeneracy of a one- 
dimensional class D topological superconductor, distin- 
guished by the global fermion parity. Using the TR trans- 
formation of 9, one can see that \n) \2m — n) under 
TR. Therefore except the states |0) and |m), all other 
states are not TR invariant. 

One may wonder whether the 2m-fold degeneracy is 
protected, and if that is the case, what physical property 
protects such degeneracy. First of all, the existence of 
gapless edge modes requires TR symmetry in the bulk 
of the FTF^l. Given the TR symmetry, because the FTI 
state is fully gapped, at energies well below the bulk gap, 
QP tunnelings through the bulk arc highly suppressed. 
Due to the fractional statistics of the QP in FTI with sta- 
tistical angle ±7r/m, any external physical perturbations 
must change the QP number of each species by m since 
they have to be local with respect to electrons. This fact, 
together with the conservation of the total fermion parity 
in the gapped superconducting systems, implies an emer- 
gent conservation of QP numbers module 2m on the edge. 
Therefore, the degeneracy of the 2m ground states with 
different QP numbers mod 2m can not be lifted by any 
physical perturbations. We therefore conclude that the 
2TO-fold degeneracy is protected by the topological order 
in the bulk FTI and the fermion parity conservation. 



Here rii is the canonical conjugate operator to the phase 
variable 9i: [9i,nj\ = i6ij. It can be regarded as the 
discrete version of the QP density. In the limit of large 
A, 9 can only takes value in a discrete set given by 0„, n = 
0, 1, ... , 2m — 1. We introduce a set of basis on each 
site and the first term cos{9i — 9i+i) becomes U^U., 
where the C/^'s have matrix representation 



[/, = diag(l,t^,w^...,t^^'"-^) 



(12) 



which is simply the operator e'*' projected onto the low- 
energy sector. 

The X]j"f term causes transitions between the \9) 
states. We rewrite this term in the \9) basis as (on one 
site) 



E 



(13) 



Here the operator V is defined as 



/O 1 
1 





Vi 







1 

0/ 



(14) 



We now make the approximation to keep only the j = 1 
term in the sum: nf oc Vi + . This is justified in the 
limit of small fluctuations of 0. It is worth mentioning 
that U and V form a 2TO-dimensional representations of 
the Weyl group algebra: 



jj2m ^ y2m ^l^yjj^ ^Jjy^ 



(15) 



So far we have succeeded in mapping the model (11) to 
the 2TO-states quantum Potts modeP^t 
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Here A cx 5^. 

We now briefly discuss the symmetry of the Potts 
model. The Hamiltonian (16 1 apparently has a global 
^2m symmetry given by 



(17) 



To see the physical meaning of the 1j2m symmetry, it is 
useful to go to a "dual" basis \n) on each site which are 
the eigenstates of n. In this basis, the global symmetry 
takes a very transparent form: Q — which is noth- 

ing but the mod 2m QP number we have mentioned in 
the previous section. 

A crucial property of the Potts model is the self dual- 
it^p^, revealed by the duality transformatiorP^J 



uL.u,. = v\v,, = u u, 



j+i> 



(18) 



which maps i?potts into itself with g ^ up to an 
overal rescaling. We can use the results in to 
obtain the explicit form of the zero modes at the special 
solvable point g = 0- Notice the product of the original 
variable and the neighbouring dual variable 7^ = Uiijj_^ 
satisfies the following algebra: 



2m 



(19) 



For m = 1 this is just the algebra of Majorana fermions. 
We call the zero modes 7's which satisfy ( 19 ) as 



zero modes. Using this result, the original model can be 
rewritten in terms of the 7^ operatorg^S 



Potts 



e'"7i,i+i72,» + Xe^"^Yi tl2A + ti-c. (20) 



The explicit expression of the 7's in terms of the original 
variables are given below: 



71, 



(21) 



which generalizes the Jordan- Wigner transformation for 
fermions. Notice that we have multiplied a phase factor 
e55s to make sure that 'yf™ — 1. 

For a finite chain with sites starting at i — and end- 
ing at i — N , = 7i,i=o and = ^2,i=N are decoupled 
from the Hamiltonian at A = 0, similar to Kitaev's Ma- 
jorana chain modeP^. One can also check that 7l,_r com- 
mute with all other terms in the Hamiltonian. In terms 
of the original bosons, the two zero modes are expressed 
as 



IL 



(22) 



When A 7^ 0, the model in general can not be solved 
exactly . However, since when A = the bulk is fully 
gapped, one can imagine slowly turning on a small but 
finite A without closing the bulk gap. The finite A phase 
should be adiabatically connected to the gapped phase 
of the FTI edge theory. During this adiabatic evolution. 



the zero modes are renormalized but remains localized 
as long as the bulk gap is not closed. In particular, we 
expect the relation (22) holds true at least in the low- 
energy sector. We also identify from ( 22 ) 



(23) 



Again we expect it to be valid universally in the gapped 
phase. 

We derive the expression of the zero modes 7^ within 
the low-energy effective Hamiltonian which takes a non- 
local form with a string operator attached. One may 
wonder whether they are local objects in terms of the 
QP operators of the FTI edge theory. Let us remark 
that the local fields Ui and Vi in the effective Hamilto- 
nian are themselves non-local in the original edge theory. 
Although a formal proof is lacking at this stage, we be- 
lieve it is highly likely the zero modes are local in terms 
of the quasiparticle operators. 



V. BRAIDING OF THE Zz^ ZERO MODES 

We now discuss the braiding statistics of the Z2,n zero 
modes. We need to specify what braiding means for the 
zero modes localized at the domain walls in one dimen- 
sion. In a closely^ related context, it has been recently 
demonstratec?^'^^ that Majorana zero modes in one- 
dimensional wires can be adiabatically exchanged either 
by forming networks out of the wires and moving the Ma- 
jorana zero modes using the "T" -junction geometry, or 
by a series of QP tunneling processes. The non-Abelian 
braiding statistics resulting from these operations takes 
the same form as the one of the Majorana fermions in 
two-dimensional -\- ipy superconductors up to an un- 
determined overall phase. This paves the way to perform 
topological quantum computing in quantum wires. There 
is no fundamental difficuties in exploiting these ideas to 
the situations at hand and therefore we can discuss the 
braiding of zero modes. 

In the following discussion, we derive possible forms 
of the braiding matrix based on general principles with- 
out appealing to details of the implementation of braid- 
ings. We will extensively use the language of the al- 
gebraic theory of anyons (namely, modular tensor cate- 
gory theory) and assume that the braiding matrices of 
the zero modes satisfy some of the basic relations origi- 
nally developed for anyons in topological phases in two 
dimensional. Whether this description applies to quasi- 
one-dimensional objects as those considered here is not 
justified a priori. It is an important question to clar- 
ify the relationship between the anyon model at hand 
and the modular tensor category theory describing two- 
dimensional topological phases. The information we need 
is the fusion algebra of the anyons as well as the mon- 
odromy equation relating the braiding matrices to the 
topological spins of the particles. 

We start from the fusion algebra of the zero 
modes. There are 2m -I- 1 types of particles, denoted 
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by (T, V'fc where k = 0,1,..., 2m — 1. Here = 1 is 
the vacuum. The nontrivial part of the fusion algebra is 
given by 



Cr X O- = V^fc, 0- X = 



k=0 



(24) 



iPp X Ipq = 1pp+q mod 2m- 



Here a represents the zero mode and ipn is the state 
with n QPs, which are Abehan with quantum dimension 
1. The vaccum fuses triviaUy with everything else. The 
fusion of two ct's follows straightforwardly from the topo- 
logical ground state degeneracy discussed in Sec. HI. The 
fusion of ip particles follows from the fact that they are 
labeled by the number of QP's mod 2m in the super- 
conducting segment of the FTI edge, and thus naturally 
form a structure. 

Using the language of modular tensor category, we 
have the following basis data from the fusion algebra: 



N±=dimVt-=l. 



(25) 



We have already counted the quantum dimensions of the 
anyons: da = V2to, d^^ = 1. 

Now let us consider braiding a pair of a particles 
labeled as 1 and 2 and for simplicity we assume that 
physically they belong to the same topological segment. 
As discussed in the previous analysis, such a segment 
has in fact only one ground state within a given mod 
2m QP number sector. This translates into the one- 



dimensionality of the fusion space as given in ( 25 1 . There- 



fore the result of an adiabatic exchange is essentially 
Abelian: the ground states in different mod 2m QP num- 
ber sector can not mix and each acquires an Abelian 
Berry phase. Let us assume that the ground state with 
mod 2to QP number n acquires a phase a„. In more ab- 
stract terms, the i?-matrix encoding the effects of braid- 
ing i?^'^ is actually one-dimensional since Nf^ = 1. 
Choosing a proper normalization, it is exactly the phase 
factor e*"". 

Without knowing how the braidings are implemented, 
one can actually determine the i?-matrix to a large extent 
using the monodromy equation in our case when fusion 
spaces are all one-dimensional. The operation of two con- 
secutive braidings, being equivalent to moving one anyon 
around another, is termed as a monodromy. It is known 
that the monodromy is fully characterized by the topo- 
logical spins of the anyons. Using the fusion algebra (24 1, 



we have the following relation between the i?-matrix and 
the topological spins of the anyons a and V'n^ 



(26) 



To figure out the topological spin of ipm we notice that 
can be regarded as a composite of n QP. Each QP, 
being an Abelian anyon in a = l/w Laughlin FQH 
liquid, has topological spin e~ where s G Z. Thus the 



topological spin 9^^ = e "« Since (i?,",," )' = e 
we find 



aa\2 ^2ia^ 



(27) 



Therefore, we have determined e*"" up to a global phase 
±0~^. This result is derived from very general consider- 
ation and should be independent of the implementations 
of the braiding operations. One can see that if m = 1 
and s = 1, (27) gives the well-known Berry phases for 



the even and odd parity ground states of two Majorana 
zero modeJ^ if we choose all the unspecified sign factors 
in ([27]) to be 1. 

We go on to consider four zero modes 71,2,3,4 arranged 
such that -fijj = e^jjji for i < j. We assume that 
71 and 72 belong to the same topological segment and 
the same to 73 and 74. Naturally we choose the Fock 
basis formed by the degenerate ground states of the two 
segments, labeled by mod 2m QP number Qi = Q12 oc 

7l72,Q2 = (934 OC 7^74: 1^1,^2) = 7"l73"|0), 711,712 = 

0, 1, ... , 2m — 1. Braiding 71 and 72 (or 73 and 74) just 
generates a Abelian phase on the Fock states. The non- 
trivial one is braiding 72 and 73 which belong to different 
topological segments. To derive the non-Abelian Berry 
phase acting on the ground state manifold, we make use 
of the Z2m conservation of QP number. After braiding, 
the QP number Q'l cx 7273 as well as Q'2 oc 7^74 should 
remain the same. This implies that if one rotates the 
Fock basis to one in which Q'^ 2 becomes diagonal, the 
braiding results in an Abelian Berry phase as discussed 
above. Then one can unwind the basis transformation to 
obtain the braiding matrix in the original basis. Simi- 
lar procedure can be carried out for any 27i zero modes, 
although our method for determining the braid matrix 
quickly becomes very inefficient. 

So far our discussion of braiding has been general. A 
realistic scheme could possibly be developed exploiting 
the idea of measurement-only topological quantum com- 
putatiorPi^, in which braidings are carried out by a se- 
ries of measurements of the fusion outcomes for pairs of 
anyons in the presence of additional topological qubits 
that are properly initialized. 



VI. FRACTIONAL JOSEPHSON EFFECT 

In this section we discuss the manifestation of the topo- 
logical degeneracy in Josephson transport. We show that 
the Josephson current has a period of Anm, thus gener- 
alizin g the An Josephson effect in topological supercon- 
ductoil™. 

Before we proceed, it is useful to understand the 
gauge transformation on the state \n). Assuming that 
we change the phase x by 27r: x ~^ X + 27r, then 



Therefore 



implying 



^3 ^ "J ^ m 
7 — > OJJ. 

In terms of the localized modes, we can write down the 
general form of the effective Hamiltonian in the tunneling 
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regime to leading order in the tunneling amplitude t: 

i^cff = r(x)7l7fl + h.c.. (28) 

Here 7^ and jji are localized zero modes at the two ends 
of the junction. Notice that here it is crucially important 
that the Josephson current is carried by quasiparticles 
instead of electrons. Although the junction is gapped out 
by TR-breaking field (e.g. ferromagnet), it is still part of 
the fractional quantum spin Hall fluid. Since they carry 
charge e/m, the period of the Josephson current should 
become 4Trm. In the following we will provide a more 
rigorous calculation of the ATrm period of the Josephson 
current. 

The DC Josephson current is then given by 



Hx) = 2e 



d(geff) 

dx 



4e5R[;|^(7l7«)]. 



(29) 



Here {j^'^r) is the conserved mod 2m QP number. 
Therefore the periodicity of / is completely determined 
by the periodicity of the phase-dependent coupling r(x). 

Now we consider the coupling r(x). We assume a 
gauge choice in which superconducting phase to the left 
of the junction is fixed and the the phase on the right 
is X- Now we increase the superconducting phase by 27r: 
X ~> X + 27r. As we just demonstrate, the zero modes 
7l — ^ 7l, 7/f — > ujjR. Equivalently, this extra phase fac- 
tor can be absorbed into the coupling T: 



r(x + 27r) =c^r(x). 
Consequently, the Josephson current 

/(x + 2^)-4e3?Lr'(x)(7l7fl) 



(30) 



(31) 



which is in general different from /(x)- It is easy to 
see that for /(x) to return to its initial value, x has to 
advance by Airm since w^™ = 1. So in summary we have 
derived that 



/(X + 27r) ^ /(x), Hx + 47rm) = /(x). 



(32) 



Although we have not determined the precise functional 
dependence of I{x), it is sufficient to conclude that the 
Josephson current has a period ATrm. 

We now briefly discuss the AC Josephson effect. The 
conventional AC Josephson current at finite voltage V 
has frequency '^j^- In our case, since the Josephson cur- 
rent is transported by a single fractionalized QP, the 
AC Josephson current has frequency When the 

Josephson junction is irradiated by microwaves with fre- 
quency w, the Shapiro steps in DC current are observed 



at Vn = 



nvahui 



i.e. we only see Shapiro steps at multi- 



ples of 2m. However, the missing steps are filled in by 
higher-order tunneling terms in the junction. A recent 
proposal of a three-leg Shapiro-step measurement may 
also be adopted here to overcome this obstacl^^. 



VII. CONCLUSION AND DISCUSSION 

In conclusion, we have investigated the exotic gapped 
phase formed on the edge of FTI under superconducting 
proximity effect. This phase is characterized by topolog- 
ical degeneracy determined by the electron fractionaliza- 
tion. The degenerate ground states are distinguished by 
the mod 2m QP number. We relate the degeneracy to 
zero modes on the boundary of the fractionalized 
topological superconducting region. We derive an effec- 
tive Hamiltonian to describe low-energy physics which 
maps to a quantum Potts model and allows for explicit 
construction of the zero modes. We then discuss the 
braiding of these zero modes. We also propose fractional 
Josephson effect as a signature of this unusual phase. 

It will be very interesting to work out the complete al- 
gebraic description of the non-Abelian excitations found 
in this work, in the framework of tensor category theory. 
One may wonder, based on the relation to the chiral Potts 
model, that the zero modes discussed in this work can be 
regarded as one-dimensional analog of parafermions 
arising in, e.g. Read-Rezayi fractional quantum Hall 
states^S. However, as non-Abelian anyons, they have 
apparently different fusion rules and braiding statistics 
and therefore should not be confused with each other. It 
is indeed true that both of them have close relation to 
Z2„i clock models. The wavefunctions of parafermionic 
fractional quantum Hall states are constructed from the 
correlation functions in parafermion conformal field the- 
orj03 which describes the critical point of certain 2,2m- 
symmetric statistical model. The zero modes studied in 
this work are related to the ordered phase of chiral Potts 
model. In both cases, the zero modes and parafermions 
are expressed by the order and disorder variables of a 
Z2m-symmetric modeP^. However, they should be placed 
properly in very different phases of the model. 

We also notice that non-Abelian anyons with quantum 
dimension d = \/N where > 1,N € Z have been found 
as lattice dislocations in Z^v gauge theorjHH and also 
lattice fractional quantum Hall state^^. The topologi- 
cal field theory of non-Abelian phases containing these 
anyons has not been found. Some candidates include the 
§0(A^)2 Chern-Simons theory^o and a U(l) x U(l) x Z2 
Chern-Simons theory studied by Barkeshli and WerP^ 
may be relevant in this context. 

We now discuss possible future directions. We have 
considered the edge states of 2D FTI. A natural ques- 
tion is whether the approach taken in this work can be 
generalize d to fractional topological insulators in three 
dimension^?'-? the surface states of which are fraction- 
alized, helical non-Fermi liquid. In contrast to the 2D 
case where at least theoretically the existence of FTI 
has been firmly established, in 3D the situation becomes 
much more complicated and so far apart from one ex- 
actly solvable modeP^, other attempts are all based on 
parton constructions and gauge theories. However, one 
could expect that superconducting proximity effect on 
such exotic surface states would lead to very interesting 
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phenomena. In particular, one can suspect that there 
may be unusual zero modes localized in the cores of su- 
perconducting hc/2e vortices. 

Note added. Upon finishing the manuscript we become 
aware of several preprints "^"^ on closely related topics. 
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